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Abstract 

We investigate the recently proposed nonlinear equation for the unintegrated gluon 
distribution function which includes the subleading effects at small x. We obtained 
numerically the solution to this equation in (x, k) space, and also the integrated gluon 
density. The subleading effects affect strongly the normalization and the x and k 
dependence of the gluon distribution. We show that the saturation scale Qs{x), which 
is obtained from this model, is consistent with the one used in the saturation model 
by Golec-Biernat and Wiisthoff. We also estimate the nonlinear effects by looking at 
the relative normalization of the solutions to the linear and nonlinear equations. It 
turns out that the differences are quite large even in the nominally dilute regime, that 
is when ^ Q^. Finally, we calculate the dipolc-nuclcon cross section. 
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1 Introduction 



The knowledge of the QCD dynamics at high energies is essential in understanding the 
hadronic interactions studied at current (HERA, Tevatron) and future (LHC) accelerators. 
Parton distributions extracted from HERA ep collider will be used in the description of the 
hadronic processes studied at LHC. It is very important to know these parton distributions 
with very high accuracy and, perhaps what even more important, to be able to estimate 
possible uncertainties which may emerge when extrapolating to the kinematic regime of 
LHC. A lot of effort is currently devoted to extracting the parton distributions with very 
high precision ^ |2j from the available experimental data. 

In principle two different frameworks can be used for calculating the parton distribu- 
tions. The standard one is based on the DGLAP evolution and collinear factorization. In 
the high energy limit it is also possible to use the fc-r factorization 3 in which the QCD 
interaction is described in terms of the quantity which depends on the transverse momen- 
tum of the gluon i.e. the unintegrated gluon distribution. An equation which governs the 
evolution of this distribution is the BFKL equation [1] . Its well known solution leads to a 
very strong power growth of the gluon density with energy: ~ s"^ where A = 4:\D.2asNc/TT 
is the BFKL intercept in the leading logarithmic approximation in powers of Ug In 1 /x 
(LLX). Next to leading order corrections to BFKL decrease the rate of growth but do 
not change the power behavior of the gluon distribution. Thus the growth of the resulting 
hadronic cross section has to be eventually tamed in order to satisfy the unitarity bound 

i- 

The perturbative parton saturation, first discussed in a pioneering paper [Jj, is a phe- 
nomenon which slows down the rapid growth of the partonic densities. It is believed that 
it leads to the restoration of the unitarity of the scattering matrix^ . When the density of 
gluons becomes very high, the gluon recombination has also to be taken into account. This 
leads to a modification of the evolution equations and their solution results in the satura- 
tion of the gluon density. In the high energy limit, the parton saturation is described as 
an infinite hierarchy of the coupled evolution equations for the correlators of Wilson lines 
[0]. It is equivalent to the JIMWLK functional equation derived within the theory of 
the Color Glass Condensate (11] . In the absence of correlations, the first equation in the 
Balitsky hierarchy decouples, and is then equivalent to the equation derived independently 
by Kovchegov 12 within the dipole formalism jl3j . 

It is desirable to have a formalism which embodies the resummation of the subleading 
corrections in Inl/x and still contains the saturation effects. Various attempts in this 
direction already exist, see for example JH-|13|- In this contribution we analyse in better 
detail the nonlinear equation for the unintegrated gluon distribution function proposed in 
|19l l2Uj . Its linear term is formulated within the unified BFKL/DGLAP framework 
and the nonlinear term is taken from the Balitsky-Kovchegov equation. We find that the 

^Unitarity or Froissart bound is valid with respect to the whole QCD, whereas parton saturation is a 
perturbative mechanism. As discussed in apart from the saturation, the confinement is also needed to 
satisfy Froissart bound in QCD. 
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subleading corrections play an important role in the calculation of the unintegrated gluon 
density since they reduce the value of the intercept and lower the overall normalization 
of the solution. We also study the differences of the solutions in the linear and nonlinear 
case. It is interesting that these differences become amplified in the case of the integrated 
gluon density xg{x,Q'^). The behaviour of the saturation scale is controlled by the value 
close to the intercept of the solution of the linear equation. In our case this value is equal 
to the one suggested by HERA data, Qs ~ exp(Ay) with A ~ 0.3. However, a more 
detailed analysis shows that even though saturation scale seems to be rather low, the 
actual numerical differences between the linear and nonlinear solutions are much bigger. 
In other words, the effect of nonlinearity on the overall normalization of the solution can 
be present already at the scales exceeding the saturation scale. 

The outline of the paper goes as follows: in the next Section we introduce the Balitsky- 
Kovchegov equation and the formalism which enables to write it in terms of the un- 
integrated gluon distribution. In Sec. El we the recall basic ingredients of the unified 
BFKL/DGLAP framework and, following |191 OU] we formulate the modified Balitsky- 
Kovchegov equation. 

In Sec. m we perform numerical analysis of this equation. We present its solution, 
i.e. the unintegrated gluon distribution as well as the integrated gluon density. Then, 
we perform the analysis of the saturation scale Qsix) and try to quantify the importance 
of the saturation effects by looking at the difference between the linear and nonlinear 
solution. In Sec. 5 we present the results for the dipole cross section a{x,r) and compare 
them to the Golec-Biernat and Wiisthoff 1^ parametrisation. We summarize our study 
in the last section. 



2 The Balitsky-Kovchegov equation 

In the dipole picture j^H] one can view the deep inelastic scattering process as a formation 
of the qq dipole, followed by the scattering of this dipole on the target. In the high energy, 
s ^ ^ ^QCD' regime these two processes are factorized, and the total 7*A^ cross 
section^ can be written as 

a^J(x,Q2) = J cfhcfirdz\^T,U^,Q^,z)\^ N{r,h,x) , (1) 

where —Q^ = q^ is the photon virtuality squared, and x ~ Q"^ / s is the usual Bjorken 
variable. The quantity ^t,l{^-,Q'^ ■, z) is the photon wave function which depends on the 
virtuality and the size of the dipole r, as well as the longitudinal fraction z of the 
photon momenta carried by the quark. Subscripts T and L denote the transverse and 
longitudinal polarisation of the incoming photon, respectively. iV(r, b, x) is the amplitude 
for the scattering of the dipole at impact parameter b on the target. It contains all the 
information about the interaction of the dipole with the target. 

''A'' being a target, nucleon or nucleus. 
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The Balitsky-Kovchegov (BK) equation |31 112j is the non-Unear equation for the am- 
phtude A^(r, b, x) 



dN{r,h,x) 
dlnl/x 



= a 



■s 



I 




r + r 



2 




- A^(r,b,x) -N{r',h + 



r + r 



2 



where = ""J^" . The Unear term on the right-hand-side of (jSJ is equivalent to the BFKL 
equation in the coordinate space, whereas the nonhnear term is responsible for the gluon 
recombination. It has been shown |4Uj to be equivalent to the triple Pomeron vertex 
[41 j . This equation has been independently derived in the dipole picture and from 
the Wilson's operator expansion [5]. In the latter case, the equation © is just the first 
member of the infinite hierarchy which decouples in the absence of correlations. 

In Eq. ((21) there is a nontrivial interplay between the two sizes: the dipole size r 
and the impact parameter b. The exact solution, recently studied in j27j (and in j28j 
with a modified kernel) has a very complicated b and r dependence which comes as 
a consequence of the conformal symmetry of this equation. Solutions to this equation 
simplified ignoring the impact parameter dependence have been extensively studied both 
analytically UHl OOl 1^ and numerically EH E31 OH • Here we are interested in 

the unintegrated gluon density f{x,k'^) averaged over the impact parameter b . Following 
jl9l l2Uj we make an ansatz that this dependence factorizes 



where R is the target size in the impact parameter. 

We are fully aware that the assumption is crude, since it implies an approximation 
of an infinite and uniform target. To obtain the full b dependence one should consider 
the exact equation 

One can now transform the equation @ in the momentum space, 



iV(r,b,x) 



n(r, x) S{h) 



(3) 



with the normalization conditions on a profile S{h) 




(4) 




(5) 



and taking 



$(l,b,j;) = ^{l,x)S{b) . 



(6) 
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We neglect the angular dependence and assume that the functions n and (p depend only 
on the absolute values of r = |r| and / = |1|. The relations (|3I5I6() enable to write the 
equation in the following form jl^ 



d(f){l,x) _ 
dlnl/x ^ 



K^(l)-^^\l,x) 



(7) 



where we have integrated both sides of over cPh. Note that while A^(r, b, x) and 
^{l,h,x) are dimensionless, the functions n{r,x) and (j){l,x) have dimension [-^^ — j] due 
to the definitions Q and ©• The operator K is the BFKL kernel [1] in momentum space 
in the LLx approximation. 

Let us now explicitly show how to find the relation between (p{l, x) and the unintegrated 
gluon distribution f{x,k'^) defined through 

xgix,Q') = J ^fix,k'), (8) 

with xg{x, Q^) being the integrated gluon density. The unintegrated gluon distribution is 
related to the dipole cross section 

a{r,x) = —J-[l- Mkr)]asf{x,e) , (9) 

which in turn can be obtained from the amplitude A^(r, b, x) by performing the integration 
over b 

a{r,x) =2 j <fhN{r,h,x) . (10) 
Using ©,0 and one obtains 

^^^^""^ = 1 1 I ^[^-Mkr)]aJ{x,k'). (11) 

Integrating over angles yields 

<P{l,x) = — ^ -Mlr)[l - Mkr)]asf{x,k^) , (12) 
and the integral over r gives 

7r2 dk^ /k^\ „, ,2, 



cP{l,x) = —J^^ —lnl^^]asfix,e). (13) 



Now we need to invert the operator 



(where g{k'^) is a test function) to get the expression for f(x,k). Multiplying both sides 
of (12) by and performing the Mellin transform with respect to we obtain 
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(/>(7,x) = / din^(l3{l,x){V) 



.2 roo ^1.2 



"•"'/Wtt^. (15) 



Nc ""(1 + 1 
and equivalently 

/(7) = ^(7 + l)V(7,^). (16) 
The inverse Mellin transform gives 

/(-,/^) = A / ^{ir''{i+im%^) = ^,{i-i'^,?m,x). (17) 



This relation between functions / and (p has been first derived in [22] and also on in 

[111201. 

3 The non-linear equation for the unintegrated density 

The relation ()17() allows us to transform Q into an equation for the unintegrated gluon 
distribution 

df{x,k^) _ a^N,^^^ r dk'^ f f{x,k'')-f{x,e) ^ f{x,e) \ 



dlnl/x vr Jkf, k'^ [ \k'^ - k'^\ [Ak'^ + k^f2 

2 d k 



■"■^''-^ dk^J m 



2 r 2 r ^nr, 7T /9 / T /9 \ T 2 



It is BFKL equation supplemented by the negative nonlinear term. 



(18) 



3.1 A partial resummation of the NLLx corrections 

Equation (|18j) contains the BFKL kernel at leading logarithmic (LLx) accuracy. This is 
a coarse approximation as far as a description of the HERA data is concerned. It is well 
known ^ that the NLLx corrections to the BFKL equation are quite large. To make the 
equation more realistic, it was proposed jl9l l2Uj to implement in the linear term of 1)18^ a 
unified BFKL-DGLAP framework developed in [2^. In this scheme |23], the BFKL kernel 
becomes modified by the consistency constraint [HSl ESI 

k'^ < k^/z, (19) 

imposed onto the real-emission part of the kernel in Eq. H18() 

r /(|,fc^2)e(^-fc^2) _ _^(|^fc2) ^ /(f,fc^) I 

/c'2 \ \k'^-k^\ |4A;'4 + A;4|I J ' 
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The consistency constraint H19() resums a large part of the subleading corrections coming 
from a choice of scales in the BFKL kernel [371 138j . Additionally, the non-singular (in x) 
part of the leading order (LO) DGLAP splitting function is included into the evolution 



J X 'J X "/ 'J X 

2N 

P99i^) = PyA^) - -f ■ (22) 

Additionally, we assume that in our evolution equation ag runs with scale A;^ which is yet 
another source of important NLLx corrections. The final improved nonlinear equation for 
the unintegrated gluon density is as follows 

/(x,fc2) ^ /(0)(^,A;2)+ 

^ a^(m^2 r /■ dk'-' f /(f,fc^^)e(f -fc-^) - /(f,fc^) ^ /(f,fc^) \ ^ 
TT L z kl k'"^ I \k'^-k^\ |4A;'4 + A;4|I J 



where 



J X J 



TT 

d^ i^.L ~ 



asik'^)lnr^]f{z,k'^) 



fc2 



k"^ ' VA:2 



(23) 



In [25 the inhomogeneous term was defined in terms of the integrated gluon distribution 



f"^ix,k') = dzP,,{z)-g [^-^X) (24) 

taken at scale k"^ = IGeV^. This scale was also used as a cutoff in the linear version of 
the evolution equation ()23() . In the linear case this provided a very good description of 
F2 data with a minimal number of physically motivated parameters, see j24| . The initial 
integrated density at scale k^ was parametrised as 

xg{x,kl)=N{l-xY , (25) 

where = 1.57 and p = 2.5. 

Let us finally note that in this model only the linear part of the BK equation has 
subleading corrections. We do not know yet how to include these corrections in the 
nonlinear term. This would require the exact knowledge of the triple Pomeron vertex [IJ 
at NLLx accuracy, which is yet unknown beyond the LLx approximation. 



4 Numerical analysis 

4.1 The unintegrated and integrated gluon density 

In this section we recall the method of solving Eq. H23|) and we present the numerical 
results for the unintegrated gluon distribution function f{x,k'^) and the integrated gluon 
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density xg{x,Q'^). The method of solving ^l'd\ . developed in [20], relies on reducing it to 
an effective evolution equation in In 1 /x with the boundary condition at some moderately 
small value of x (i.e. x = xq ~ 0.01). 

To be specific, we make the following approximations: 

1. The consistency constraint j z — k''^) in the BFKL kernel is replaced by the 
following effective {z independent) term 



e{e/z - k'^) ^ e{e - y') + ( 



12 \ 



\k'^ 



e( 



u'2 



(26) 



This is motivated by the structure of the consistency constraint in the moment space, 
i.e. 



'2 k^) 



2. The splitting function is approximated in the following way 

'-^ dz 



^ [zPg,{z) - 2N,]f [-^,k'') ^ Pgg{u: = 0)fix,k") 



where Pgg{uj) is a moment function 



'^z-[zPgg{z)-2N,] 



and 



Pggi^ = 0) 

nd 

of Pgg{uj) around u = 0, see P^j . 

Using these approximations in H23|l we obtain 

df{x,k^) ^ 
91n(l/x) 



(27) 

(28) 

(29) 
(30) 
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ggyUJ = [)) = -—. 

This approximation corresponds to retaining only the leading term in the expansion 



a. 



2 A;'2 





2)iv. f dk'^ r / ^'^) [0 (^^ - k") + (1^)"' e (P - k'-)\ - f (x, k-) 



\k'^-k^\ 



fix,k^ 



|4fc'4 + k^ 



1 - k' 



IT 



d V k^ 



dk'^ 



2 A;'2 





dk"^ i?2 



oo J1J2 



k2 



dk' 
1^ 



■as{k'')ln 



f{x,k'^ 
fk'^ 



\k^ 



fix,k 



I2\ 



■ (31) 



First, the equation ()31() was solved with the non-linear term neglected starting from 
the initial conditions at x = 10"^ given by (|24j) . The parameter was adjusted in such 
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Figure 1: The unintegrated gluon distribution f{x, k'^) obtained from Eg. it^j|) as a function 
of X for different values fc^ = 10^ GeV^ and P = 10^ GeV^ (left) and for k'^ = 5 GeV^ 
and fc^ = 50 GeV^ (right) . The solid lines correspond to the solution of the nonlinear 
equation i^Jl]) whereas the dashed lines correspond to the linear BFKL/DGLAP term in 



a way that the solution of the linear part of matched the solution of the original 
equation in the BFKL/DGLAP framework "24 . This procedure gives lo^q = 0.2 and the 
solution of the linear part of (I.Slj) reproduces the original results of [23] within 3% accuracy 
in the region 10"^ > x > 10~® and 2 GeV^ < k"^ < 10^ GeV^. This matching procedure 
has also the advantage that the quark contribution present in the original BFKL/DGLAP 
framework is effectively included by fitting the value of Weff- The full non-linear equation 
pif) was then solved using the same initial conditions and setting R = A GeV~^. 

In Fig. n we plot the unintegrated gluon distribution function as a function of x for 
different values of k'^. This figure compares the results of two calculations, based on 
the linear and nonlinear equations. The differences are not large, however there is some 
suppression due to the nonlinearity at smallest values of x < 10~^. 

The subleading corrections strongly decrease the value of the intercept with respect 
to the LLx value and the nonlinear term becomes important only at very low values of 
x. As is evident from Fig. [21 the subleading corrections cause a large suppression in the 
normalisation, also at moderate values of x. This is due to the fact that the non-singular in 
X part of the Pgg splitting function was included into the evolution. This term is negative 
and is important at large and moderate values of x. 

The same conclusions can be reached by investigating the plots in Fig. |21 where the 
unintegrated density is shown as a function of the transverse momentum fc^ for fixed 
values of x. The nonlinear effects seem to have a moderate impact in that region. On 
the other hand the subleading corrections are substantial. For example, at x = 10~^ and 
k'^ = 10 GeV^ the reduction in magnitude of the unintegrated gluon density is about 25%. 

In Fig. |1] we show the integrated gluon density given by Eq. © . The change from the 
power behaviour at small x is clearly visible in the nonlinear case. Also the differences 
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Figure 2: The same as Fig. ^ but now the modifed BK equation H31|) (soUd lines) is 
compared with the original BK equation (jlBj) without subleading corrections (dashed 
lines) . 



between the distributions in the linear and nonlinear case seem to be more pronounced for 
the quantity xg{x,Q'^). This is due to the fact that in order to obtain the gluon density 
xg{x,Q'^) one needs to integrate over scales up to including small values of k'^, where 
the suppression due to the nonlinear term is bigger. 



4.2 The saturation scale Qs{x) 

In order to quantify the strength of the nonlinear term, one introduces the saturation 
scale Qsix). It divides the space in (x,/c^) into regions of the dilute and dense partonic 
system. In the case when k'^ < Q'l{x) the solution of the nonlinear BK equation exhibits 
the geometric scaling. This means that it is dependent only on one variable N(r, x) = 
N{rQs{x)), or in momentum space 4'{k,x) = (t){k / Q s{x)) . Recently, an analysis of the 
saturation scale in the case of the model with resummed NLL BFKL has been performed 
jl8j . There, the saturation scale was calculated from the formula 

d'ye 1 - 7c ' 

which has been first derived in [7] by the boundary condition of the wave front. Formula 
H32() has been later rederived in [421 l3Uj . The effective Pomeron intercept ojs is a solution 
to the equation 

ujsil) = asx{l,^s) , (33) 

where ^(7)'^) is the kernel eigenvalue of the resummed model. In our case the eigenvalue 
has the following form 

X(7, w) = 2^'(1) - vI/(7) _ ^(1 - 7 + ^) + -Pgg{oj) . (34) 

7 

The solution for the saturation scale obtained from solving H32I33|1 using eigenvalue H34|) is 
shown in Fig. and gives A = '^^f^ = 0.30,0.45,0.54 for three values of = 0.1,0.2,0.3, 
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Figure 3: The unintegrated gluon distribution f{x, k'^) as a function of k"^ for two values 
of X = 10~^ and 10~^ . Left: solid lines correspond to the solution of the nonlinear 
equation \31)) whereas dashed lines correspond to linear BFKL/DGLAP term in 
Right: solid lines correspond to the solution of the nonlinear equation h'Jl]) whereas dashed 
lines correspond to the solution of the original BK equation without the NLLx modifications 
in the linear part il^). 



respectively. These results are similar to those obtained in ^H]. We compare our results 
with the saturation scale from the Golec-Biernat and Wiisthoff model. Normalisation of 
the saturation scale is set to match GBW saturation scale at xq = 0.41 x 10~^. 

The saturation scale Qsix) can be also obtained directly from the numerical solution 
to the nonlinear equation by locating, for example, the maximum of the momentum distri- 
bution of the unintegrated gluon density in the spirit of method presented in Ref. . For 
the purpose of phenomenology we attempt here to estimate the effect of the nonlinearity 
in a different, probably more quantitative way. We study the relative difference between 
the solutions to the linear and nonlinear equations 

/i-(:E,a(x,/3)2) 

where /? is a constant of order 0.1 — 0.5. Since this definition of the saturation scale is 
different from the one used in the literature and is likely to posses different x dependence, 
we denote it as Qg- In Fig. I^I^left) we show a set Qs which are solutions to Eq. (|35|) 
for different choices of /? together with the saturation scale calculated from the original 
saturation model by Golec-Biernat and Wiisthoff |2^. Solid lines given by Eq.(|35|) show 
where the nonlinear solution for the unintegrated gluon starts to deviate from the linear 
one by 10%, 20%, . . . , 50%. It is interesting that contours Qsix) defined in (|35|) have much 
stronger x dependence than saturation scale Qsix) defined by Ea. H32() and the one from 
GBW model. In particular Qsix, (3) > Qsix) for given x (at very small values of x). 
This might be a hint that saturation corrections can become important much earlier (i.e. 
for lower energies) than it would be expected from the usual definition of the saturation 
scale Qsix). In Fig. Efright) we also show contours in the case of the integrated gluon 
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Figure 4: The integrated gluon distribution xg{x, k'^) as a function of x for values of = 
lO^GeV^ and = iq3 ^i^flj and for Q"^ = 5 GeV^ and = 59 (right) obtained from 
integrating fix. k!^)Eq. iSl]). Dashed lines correspond to solution of linear BFKL/DGLAP 
evolution equation. 

distribution function, that is the solution to 1)35(1 with /(x, k'^) replaced by xg{x., Q^). As 
already seen from the previous plot, Fig. ^ the differences in the integrated gluon are 
more pronounced. For example in the case of = 25GeV^ and x ~ 10~^ — 10~^ we 
expect about 15% to 30% difference in the normalization. Again, by looking solely at 
the position of the critical line, one would expect the nonlinear effects to be completely 
negligible in this region since at x = 10~^ the corresponding Ql{x) ~ 2.8GeV^ (taking 
Qlix) = Q1q{x/xq)~^ with normalization Ql^ = IGeV^ at xq ~ 4 x 10"^ and A ~ 0.28, 
j25j). This rough analysis shows that one cannot think of saturation scale as a definite 
and sharp border between very dilute and dense system. The transition between these 
two regimes appears to be rather smooth and the nonlinear term of the equation seems 
to have quite a large impact on the normalization even in the 'linear' regime defined as 



In practice, the estimate of the saturation effects is even more complicated since the 
unintegrated gluon density has to be convoluted with some impact factor, and the inte- 
gration over the range of scales must be performed. 

4.3 Dipole cross section a{r, x) 

It is interesting to see what is the behavior of the dipole cross section a[r,x) obtained 
from the unintegrated gluon density via Eq. ©. In this calculation we assume that as is 
running with the scale k"^. 

Calculation of the dipole cross section requires the knowledge of the unintegrated gluon 
density for all scales < /c^ < 00. Since in our formulation the unintegrated gluon density 
is known for k"^ > k^ we need to parametrise f{x,k'^) for lowest values of k"^ < k^. We use 
the matching condition 



Q^->Ql{x). 




(36) 
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Figure 5: Saturation scale obtained from Eqs. \3}^'j'J\) solid lines, compared with saturation 
scale from GB W model iZ5^ . 

and following |43| we assume that f{x,k^) ^ for low k^. This gives (compare Ea. (|25l) 
) 

f{x,k^) =4N{l-xYk^ . (37) 

In Fig. Q we present the dipole cross section as a function of the dipole size r for 
three values of x = 10^^, 10~^ and 10^^. For comparison we also present the dipole cross 
section obtained from GBW parametrisation. To be self-consistent, we cut the plot at 
r = 2 GeV~^ because we assumed in the derivation of formula (O that the dipoles are small 
in comparison to the target size (we assume proton radius to be 4 GeV~^). This cut allows 
us to obtain a model independent result since we observe that different parametrisations 
of f{x, k"^) for k"^ < fcg give essentially the same contribution for r < 2GeV~^. 

We observe that our extraction of the dipole cross section gives similar result to the 
GBW parametrisation. The small difference in the normalisation is probably due to the 
different values of Xg which probe the gluon distribution (or alternatively the dipole cross 
section). In the GBW model the dipole cross section is taken at the value Xg = x which 
is the standard Bjorken x = Q'^/2p ■ q. On the other hand, in the formalism presented in 
Ref. [3., one takes into account the exact kinematics (energy conservation) in the photon 
impact factor. It is a part of the subleading effect in the impact factor and it increases the 
value of Xg ~ 5x. Therefore, in our formalism the normalization of the unintegrated gluon 
is increased so that the convolution with the impact factor and the resulting structure 
function remains the same. 
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k^(GeV^) Q"(GeV^) 

Figure 6: Solid lines show contours of constant relative difference between solutions 
to linear and nonlinear equations, Eq. I\'^5\) . Lines from bottom to top correspond to 
10%, 20%, 30%, 40%, 50% difference. Left: Contours in the case of the unintegrated gluon 
distribution f{x,k'^); right: contours in the case of the integrated gluon distribution 
xg{x,Q'^). Dashed line in both case corresponds to the saturation scale from Golec-Biernat 
and Wusthoff model \2,^ . 



5 Conclusions 



In this paper we studied numerically the solutions to the modified BK equation in the ap- 
proximation of the infinite and uniform target. The modifications include the subleading 
corrections in In 1/x which are given by the kinematical constraint, DGLAP Pgg splitting 
function and the running of the strong coupling. Since these corrections reduce signifi- 
cantly the value of the BFKL intercept they also have large impact onto the behavior of 
the saturation scale and the normalization of the solution. For example, we find that at 
X = 10~^ and k"^ = 100 GeV'^ the normalization of the unintegrated gluon distribution 
is reduced by about 30% — 40% as compared with the solution to the unmodified BK 
equation. 

We have studied the onset of the nonlinear corrections by observing the difference in 
the normalization of the linear and non-linear solutions. We observe that even though the 
solutions are in the nominally dilute regime, the normalization of the solution to the BK 
equation can be already strongly affected by the presence of the nonlinear term. This can 
have potential impact onto the extrapolation of the parton distributions to lower values 
of X. We note however that, as long as the nonlinearities do not affect substantially the k 
and X dependence of the solution, the linear equation can probably be used with suitably 
chosen boundary conditions. 

We have also computed the dipole cross section from this model and compared it with 
the GBW. We find that both models give similar results, with some small differences which 
can be atrributed to slightly different treatment of the photon impact factor. 

Finally, we stress that although the subleading corrections are included in this formal- 
ism by using the available knowledge on the NLLx BFKL equation and the resummation 
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30^ r 




r(GeV b 

Figure 7: The dipole cross section obtained from modified BK (solid line) compared to 
GBW dipole model (dashed line). 

procedures, they are taken into account only in the hnear part of the BK equation. Con- 
sistent and complete treatment would require their inclusion in the triple Pomeron vertex 
part, which is so far known to LLx accuracy only. 
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